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Using ensembles of two, three and four spheres immersed 
in a fermionic background we evaluate the (integrated) den- 
sity of states and the Casimir energy. We thus infer that 
for sufficiently smooth objects, whose various geometric char- 
acteristic lengths are larger then the Fermi wave length one 
can use the simplest semiclassical approximation (the con- 
tribution due shortest periodic orbits only) to evaluate the 
Casimir energy. We also show that the Casimir energy for 
several objects can be represented fairly accurately as a sum 
of pairwise Casimir interactions between pairs of objects. 

PACS numbers: 03.65. Sq, 21.fO.Ma, 26.60.+C, 97.60. Jd 



In 1948 Casimir predicted the existence of a very pe- 
culiar effect, the attraction between two metallic parallel 
plates in vacuum |]l| . The existence of such an attraction 
has been confirmed experimentally with high accuracy 
only recently |^ . The origin of this attractive force can 
be traced back to the modification of the spectrum of zero 
point fluctuations of the electromagnetic field. Similar 
phenomena are expected to exists for various other (typ- 
ically bosonic) fields |^,^ and the corresponding forces 
are referred to as Casimir or fluctuating interactions. A 
related interaction arises when the space is filled with 
(noninteracting) fermions, which is particularly relevant 
to the physics of neutron stars [||J|] and quark gluon 
plasma [Q. Spin-orbit interaction will be neglected as 
well. One of the simplest cases corresponds to nuclei 
embedded in a neutron gas. These however could be re- 
placed with buckyballs immersed in an electron gas, in 
liquid mercury for example. Particularly attractive can- 
didates for the study of this type of Casimir effects in 
essentially perfect degenerate fermi systems are the di- 
lute atomic Fermi condensates 

In the case of two parallel impenetrable planes, di- 
mensional arguments suggest that the dependence of the 
Casimir energy for fermions on the distance between the 
two planes has the form 



(1) 



where fx = fi kp/2m is the chemical potential, kp is the 
Fermi wave vector and d is the distance between the two 
planes. For this simple geometry it is straightforward to 
evaluate the function F{kpd) One has to be careful 
and specify whether the calculation should be performed 
at fixed particle number or fixed chemical potential, as 
one can easily show that the Casimir energy has a differ- 
ent behavior in these two limits. 



For more complicated geometries the evaluation of 
the Casimir energy is generally a rather involved, even 
though straightforward, numerical procedure. Our main 
goal is to reach a qualitative understanding of the 
Casimir energy in the case of complicated geometries. 
We shall consider mainly two obvious limits, when the 
objects immersed in the Fermi environment are either 
much smaller or much larger than the Fermi wave length. 
The limit of small scatterers is relatively easy to treat and 
is considered mostly for the completeness of the analy- 
sis. We show here that the case of large scatterers can 
be treated quite accurately using semiclassical methods 
at practically all separations. The most important con- 
clusion we are able to draw from our study however is 
that the Casimir interaction energy in the case of more 
than two scatterers can be evaluated quite accurately as 
a sum of pairwise interactions between these scatterers. 
This conclusion comes to some extent as a surprise, since 
it is known that Casimir energy is not pairwise additive, 
in other words, the interaction among extended objects 
cannot be evaluated as a sum of pairwise interactions, 
see e.g. Ref. |,|. 

Let us consider at first the case of two impenetrable 
spheres of radius a at a distance r > 2a between their 
centers. In order to calculate the Casimir energy we 
shall represent the sufficiently smoothed fermion density 
of states (smoothing is over an energy interval larger than 
the level spacing in the volume V of the entire system): 



lis, a, r) = go{e) + gw{£, a) + gc{£, a, r). 



(2) 



where g{e, a, r) is the total fermion density of states, go{e) 
is the density of states in the absence of scatterers (ideal 
Fermi gas), gw{£,o,) is the correction to the density of 
states arising from the presence of two spheres infinitely 
apart from each other and gc{£,a,r) is the remaining 
part, which is of central interest to us here (in the follow- 
ing we shall not make explicit the a- or i — dependence, 
but show only the energy (e-) dependence). 

In the case of N scatterers the Krein formula [p|jic|] 
provides a link between the A^-body scattering matrix 
5jv(£) and the change in the density of states due to the 
presence of N scatterers, namely 



'55(e) = 5(e) - .9o(e) = tt- 
2m 
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Following Refs. |T^, the determinant of the S'Ar(£)- 
matrix can be represented as follows 
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where M(e) is a Koringa-Kohn-Rostoker (KKR) multi- 
ple scattering matrix [|l2| . determines the change 
in the density of states due to the presence of isolated 
scatterers, which in case of large scatterers is given basi- 
cally by a Weyl formula, see Kefs. 11,0 for various ex- 
amples and general formulas, gc{£)^ which is determined 
by the multiple scattering KKR-matrix M(e), vanishes 
in the limit of infinitely separated scatterers and is the 
only part of the density of states which depends on the 
relative arrangement of the scatterers. The Casimir en- 
ergy at fixed particle number can then be introduced as 
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where the omitted terms are 0{V^'^). N is the total 
number of fermions, /i and /ip are the values for the chem- 
ical potential with the scatterers at finite and infinite 
separation from each other respectively and A/c(e) is the 
relevant correction to the integrated density of states. 
Strictly speaking, the quantities g{e) and go(£) are in- 
finite, as they are proportional to the volume V of the 
entire space. This redundant divergence can be handled 
easily by considering first a very big box, the volume of 
which is subsequently taken to infinity. One can then 
show that Ec has a well defined and finite value in this 
limit. 

Using the explicit formulas for the KKR-matrix from 
Ref. |]lT| , one can compute numerically gci^) for vari- 
ous arrangements of hard spherical (or circular in 2D) 
scatterers. It is possible to obtain significantly simpler 
expressions for the (integrated) density of states in the 
limit of very small and very large scatterers. If the wave 
length A = 27r/fc (e = h^k"^ /2m) is much larger than the 
radii of the scatterers and the scatterers do not overlap 
then one can show that the KKR-matrix M{£) is given 
by (see for the analog in the 2D case) 

[M{e)U, « <5„„, - (1 - ^„„.)/„(£) '^P^'^'"'"\ (9) 

where the indices 71,^1= 1, . . . , run over the scatterers, 
r„m is the distance between the centers of the n-th and 
m-th scatterers, /n(e) is the s-wave scattering amplitude 



on the n-th scatterer. In the case of two spheres of radius 
a, with their centers r apart (r ^ a) one obtains 

Mcie) - sin[2k{r ~ a)] + O ( {ka)\ ^) , (10) 

where v is the spin degeneracy factor. The next order 
correction arises from p-wave scattering. In the case of a 
finite radius a, one can use the Gutzwillcr trace formula 
to determine this correction to the (integrated) density 
of states semiclassically {sd) jl4| 
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where the summation is over periodic orbits (po), Tppo 
and n are the period and number of repetitions of the 
primitive periodic orbit (ppo), Mpo, Spo and the 
stability matrix, classical action and Maslov index (which 
counts the number of bounces under Dirichlet boundary 
conditions) of the po [|lj]. Taking into account only the 
contribution arising from the single po of length 2(r — 2a), 
with no repetitions, one derives the following result 



A/'c(e) 



47rr(r — 2a) 



sin[2fc(r - 2a)]. 



(13) 



At large distances (r ;3> a) the leading term in both cases 
(small and large scatterers) has the same analytical struc- 
ture, apart from an overall numerical factor. The fcr ^ 1 



limit of Rel. (13) can be reproduced from the KKR- 



matrix in the case of large spheres {ka > 1) as well. One 
can expect that the semiclassical result ( p^ ) is reasonably 
accurate when the reduced action along the classical po is 
larger than unity, i.e. when Spo/h — 2k(r — 2a) > 1, and 
that this approximation should fail when the two spheres 
are very close. Surprisingly however, at smallest separa- 
tions when r — 2a <C a the semiclassical estimate is only 
about 30% lower than the exact result. For 2k{r—2a) > 1 
the semiclassical expression (13) is very close to the exact 
numerical result obtained using the KKR-matrix, if the 
spheres are large {ka > 1), see Fig.|l|. When fca ^ 1 a 
large number of partial waves contribute, which renders 
thus the semiclassical limit valid. One can show, using ar- 
guments along the lines in Ref. |l^ ] , that the contribution 
arising from creeping orbits are exponentially suppressed, 
which is intuitively expected. The contributions arising 
from repetitions of the primitive periodic orbit are rela- 
tively small, because the long orbits and the repetitions 
of a primitive orbit are less stable in 3D than in 2D. For 
this reason the simplest semiclassical approximation for 
the density of states, and consequently for the Casimir 
energy, are more accurate for spheres than for cylinders. 

Our findings suggest that for more complicated ge- 
ometries, if the curvature radii are larger than the wave 
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length, one can safely evaluate the density of states using 
the contributions arising from short primitive periodic or- 
bits only, without any repetitions. If the curvature radii 
are smaller than the wave length then the alternative 
simple approximation of small scatterers could be used. 
Using Eq. (|l3|) one can now easily derive an approximate, 
but rather accurate, expression for the Casimir energy for 
two large spheres {kpa > 1) 
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3i[2kF{r~2a)], 



(14) 



where ji(x) is the spherical Bessel function. 
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FIG. 1. The correction to the number of states Nc{e) as 
a function of r — 2a in the upper part and as a function of k 
in the lower part for u = 1 

Under the same approximations the Casimir energy of a 
large sphere at a distance r from an infinite plane reads 
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n[2kp(r-a)]. 



(15) 



One can naturally expect that if a standing wave with 
the Fermi momentum could be formed in between the 
two spheres then the total energy of the system is at 
a (local) minimum, which thus explains the oscillatory 
character of this interaction. 



FIG. 2. The relative amplitude of the contribution to 
the semiclassical density of states (Q of n-bounce periodic 
orbits for the system of three identical spheres, situated at 
the vertices of an equilateral triangle, as a function of the 
distance between two spheres r /a{r > 2a) as compared to the 
amplitude of the 2-bounce orbit, including the corresponding 
degeneracies. In the legend 123 denotes an orbit starting at 
the sphere 1, followed by a bounce off the sphere 2, then off 
the sphere 3 and ending on sphere 1, and so forth. 

Let us consider now the case of three and four spheres. 
The semiclassical picture makes particularly transparent 
the reason why, strictly speaking, the Casimir energy 
cannot be represented as a sum of pairwise interactions. 
Each primitive periodic orbit and its repetitions give rise 
to an additive contribution to the density of states, see 
Eq.(|l2|), and thus to the Casimir energy (^. For three 
or more objects there are periodic trajectories (or stand- 
ing waves) bouncing off three or more such objects and 
thus the contribution to the density of states and to the 
Casimir energy due to such orbits depends on the rela- 
tive arrangement of three or more objects, thus leading 
to genuine three and more body interactions. We de- 
termined however that the contribution of three or more 
bounce orbits to the density of states, and thus to the 
Casimir energy as well, is never dominant. (NB Our 
analysis and conclusions refer to the global domain of 
the system always and not to the fundamental domain 
or to the one-dimensional representations of a symmetry 
reduced problem.) An analysis of the stability matrix of 
an n-bounce orbit shows that its contribution to the inte- 
grated density of states at large distances is proportional 
to 1/L", where L is the length of the orbit, if all the legs 
of the orbit are comparable in length, see Fig. (|^). Any 
person who ever played pool (thus in 2D) knows instinc- 
tively that long shots arc more difficult than short ones 
and that the most difficult shots are the many-bounce 
shots. In 3D and higher dimensions orbits are typically 
more unstable than in 2D. An exact evaluation of the 
stability matrix for various periodic trajectories shows 
that even at small separations the contributions of 2- 
bounce periodic orbits dominate over those of three or 
more bounce periodic orbits. The 3-bounce orbit gives 
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the largest contribution, an approximately 10 % correc- 
tions, at r « 2.5a. As one can also judge from Figs. 

the role of the orbits bouncing among three or more 
objects is never too large. The Casimir energy for three 
identical spheres (at the vertices of an equilateral tri- 
angle) satisfies the approximate relation E3 « 3i?2 and 
correspondingly in the case of four identical spheres (at 
the vertices of a tetrahedron) one has E4 « 6E2 ~ 2£^3, 
where En stands for the Casimir energy of N spheres, 
with high accuracy if kpa ^ 1. For kpa < 1 corrections 
could reach 10% for three and 25 % for four spheres. 
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FIG. 3. The ratio of the exact Casimir energy and chem- 
ical potential Ec / fJ> for four spheres and computed as a sum 
of the exact contributions due to pairs or triplets of spheres 
for two different separations. 

We presented here results only for the symmetric ar- 
rangement of the spheres due to the lack of space. Vari- 
ous asymmetrical configurations of three and four spheres 
show the same general pattern, namely that the correc- 
tion to the integrated density of states Mc{£) can be 
represented fairly accurately as a sum of the correspond- 
ing corrections for pairs of spheres. Obviously, one can 
replace the spheres with other objects, with curvature 
radii larger than the Fermi wave length. The pairwise 
additivity of the Casimir interaction is reasonably well 
satisfied as well for the case of point scatterers, as one 



can easily check using Eq. (||). We thus conclude that 
genuine many-body Casimir interactions are relatively 
short ranged and that two-body interactions strongly 
dominate - even for small separations. 
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